Astronomy & Astrophysics manuscript no. final radiozcosmos 


© ESO 2009 


November 3, 2009 





Properties and environment of Radio Emitting Galaxies 
in the VLA-zCOSMOS survey * 

S. Bardelli 1 , E. Schinnerer 2 , V. Smolcic 3 , G. Zamorani 1 , E. Zucca 1 , M.Mignoli 1 , C.Halliday 4 , K.Kovac 5 , P.Ciliegi 1 , 
q{ \ K.Caputi 5 , A.M. Koekemoer 6 , A.Bongiorno 7 , M.Bondi 8 , M.Bolzonella 1 , D.Vergani 1 , L.Pozzetti 1 , C.M.Carollo 5 , 
T.Contini 9 , J.-P.Kneib 10 , O.Le Fevre 10 , S.Lilly 5 , V.Mainieri 11 , A.Renzini 12 , M.Scodeggio 13 , G.Coppa 1 , O.Cucciati 10 , 
S.de la Torre 14 , L.de Ravel 10 , P.Franzetti 13 , B.Garilli 13 , A.Iovino 14 , P.Kampczyk 5 , C.Knobel 5 , F.Lamareille 9 , 
J.-F.Le Borgne 9 , V.Le Brun 10 , C.Maier 5 , R.Peiio 9 , Y.Peng 5 , E.Perez-Montero 9 ' 27 , E.Ricciardelli 12 , J.D.Silverman 5 , 

M.Tanaka 11 , L.Tasca 13 , L.Tresse 10 , U.Abbas 15 , D.Bottini 13 , A.Cappi 1 , PCassata 16 , A.Cimatti 17 , L.Guzzo 14 , 
A.Leauthaud 18 , D.Maccagni 13 , C.Marinoni 19 , H.J.McCracken 20 , PMemeo 13 , B.Meneux 7 - 28 , POesch 5 , C.Porciani 21 , 
R.Scaramella 22 , PCapak 23 , D.Sanders 24 , N.Scoville 3 , Y.Taniguchi 25 , and K.Jahnke 26 



o 

CO 



(Affiliations can be found after the references) 
Received — ; accepted — 

ABSTRACT 



u 

Aims. We investigate the properties and the environment of radio sources with optical counterpart from the combined VLA-COSMOS and 
Q . zCOSMOS samples. The advantage of this sample is the availability of optical spectroscopic information, high quality redshifts, and accurate 
'_h density determination. 

' Methods. By comparing the star formation rates estimated from the optical spectral energy distribution with those based on the radio luminosity, 
| we divide the radio sources in three families, passive AGN, non-passive AGN and star forming galaxies. These families occupy specific regions of 
1 1 , the 8.0-4.5 fim infrared color- specific star formation plane, from which we extract the corresponding control samples. 

Results. Only the passive AGN have a significantly different environment distribution from their control sample. The fraction of radio-loud passive 
' AGN increases from ~ 2% in underdense regions to ~ 15% for overdensities (1+8) greater than 10. This trend is also present as a function of 
^ ' richness of the groups hosting the radio sources. Passive AGN in overdensities tend to have higher radio luminosities than those in lower density 
m : environments. Since the black hole mass distribution is similar in both environments, we speculate that, for low radio luminosities, the radio 
emission is controlled (through fuel disponibility or confinement of radio jet by local gas pressure) by the interstellar medium of the host galaxy, 
'■O while in other cases it is determined by the structure (group or cluster) in which the galaxy resides. 
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1 . Introduction intermittent mechanis m is indirectly obse rved as bubbles in the 

I " ■ intra-cluster medium (Birzan et al., 2004) or in the restar t ed ac- 

> ■ Radio sources have been traditionally classified into two broad ^ of radio sources ( Vgnturi Dallacasa & Ste fanachi, 2004; 

^ ; classes depending upon their emission mechanism: AGN-driven Bar delli et al.L l2002h. This heating emission is thought also to be 

K£ . emission has been supposed for early type galaxies, while a star ^ tQ suppress ^ staf formation and t0 be an important mec h- 

03 forming ongir^of the^dio power has been invoked for late type anifjm for creati the al color bimodality dCroton etall 

gdaxies jCondori [1223). However not all optical galaxies ex- 2QQ6) Moreovej . it is a fact that many radio AGN reside in earl 

hibit radio emission (at least at the depth of current surveys) and laxies d Led i 0W & QwenL [T99I , a class of objects which 

it is important to investigate the physical mechanism which trig- ^ ferentiall found in hi h densit environments, 
gers this emission. For AGN-induced activity, the difference be- 
tween radio quiet and radio loud galaxies is naturally connected 

with the availability of fuel for the central en gine, the black 11 has been proposed that cluster/group mergers, galaxy- 
hole mass, and its emission efficiency (see e.g. IShabala et all S alax y mergings or at least tidal disturbance caused by close 
2008, and references therein), while for star forming objects it is encounters between galaxies would increase radio activity. For 
closely related to the strength of the star formation episode. earl y ^ galaxies, these phenomena could drive gas more effi- 
Heating by a central AGN is also thought to be impor- centl y toward the center ' whlle for late ^ g a l™, these en- 
tant for the host galaxy evolution, providing the energy to stop eomrterswould trigger and/prenhance star formation (see e.g. 
the central cooling that drives gas to the central black hole in ™i ll999t | Vollmer et al! |2001|). The radio luminosity of the 
an intermittent feedback mechanism dubbed "radio mode" (see AGN also de P ends on the interaction of the jet coming from the 
iCroton etall HI |Ck)tti Ostriker & Progal HI). The same black hole wlth the surrounding gas of the hostgalaxy and/or of 
the hosting cluster (see e.g. Sha bala et al.L 120081 and references 



Send offprint requests to: S.Bardelli sandro.bardelli@oabo.inaf.it therein). For these reasons, some degree of dependence of radio 

* based on data obtained with the European Southern Observatory activity on the environment is expected. Although the problem 

Very Large Telescope, Paranal, Chile, program 175.A-0839 and VLA- is clear from the qualitative description, controversial results are 

Large Program VLA-AS0801 present in the literature and for mainly two reasons. 
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First of all, very few papers are primarly based on radio data. 
Usually, papers in the literature consider AGN selected on the 
basis of their optical ( mostly estimated from l ine diagnostics) 
or X-ray activity (see Kauffman n et al.L 12004 and references 
therein). It seems that the family of optical/X-ray AGN inhabits 
similar environments to non-active galaxies, with a preference 
for lower densities at h igher st ellar masses (Kauff mann et all 
l2004t ISilverman etail l2009bl) . However, lBestetal.1 d2005l) 
claim that the phenomena generating the emission lines and the 
radio emission are statistically independent, even if recent re- 
sults indicate that radio loud emi ssion line AGN favour higher 
densi ties than radi o quiet one s ( Kauffmann Heckman & Besi 
2008) (see also ICaputi et all |2009l for a study in the IR). 
Secondly, only specific envi ronments have been considered , 
mainly clusters or groups dMiller et all 120021: iHill & Lillvl 
I199U iGiacintucci et all 120041) . lacking therefore the coverage 
of a large dynamical range in densities. 

The first paper that studied in a complete statistical way, both 
from the radio a nd op tical side, the environment of radio galax- 
ies was that of IBestI ( 12004) . This paper use s data from 2dF 
Galaxy Redshift Survey dColless et all |2001|) and the NRAO 
VLA Sky Survey dCondon et all 1 19981) with a well defined esti- 
mator of density, i.e. the projected density of the 10 rh nearest 
neighbor. Moreover, groups (and clusters) are found with the 
standard friend-of-friends detection method. The redshift range 
is 0.02 < z < 0.1. The general result is that radio-selected star 
forming galaxies decrease in number with increasing overden- 
sity, while objects with AGN activity exhibit little dependence 
upon the local density, except at the very low densities, where 
the pr obability to find these objects decreases significantly. Best 
(2004) claims also that the larger scale (i.e. on the scales typical 
of groups and clusters) is more important than the smaller scale 
(typical of galaxy pairs or companions) to determine the AGN 
radio emission. Star forming galaxies do not e xhibit dependence 
on en vironment at scales larger than one Mpc ( Kauff mann et all 
120041) . 

In this paper we investigate the environment properties of 
the radio sources present in the zCOSMOS survey using a set of 
well defined density estimators in the redshift range [0. 1-0.9]. In 
Section 2 we present the used data and the estimators for vari- 
ous environments; in Section 3 we describe the method adopted 
to separate galaxies with radio emission induced by AGN from 
that induced by star formation and the definition of the control 
samples. In Section 4 we discuss the differences between the 
various samples, while in Section 5 we present the dependences 
on environment. In Section 6 the luminosity distributions and 
radio-optical ratio for a specific class of radio AGN are shown 
and in Section 7 we discuss the results. The adopted cosmology 
has Q,„ = 0.25,Q A = 0.75 and H = 70 km s _1 Mpc 1 . 

Throughout the paper we adopted the term "radio galaxy" 
for galaxies hosting a radio source, independently of the physical 
nature of the emission. 



2. Data and density estimators 

To select the optical sample, we used the zCOSMOS Bright 
surve y, a redshif t survey with a magnitude selection of Iab < 
22.5 (Lill y et all 120071 l2009h th at has been undertak en in the 
1.7 deg 2 of the COSMOS field dScoville et all l2007l) . imaged 
by the Hubble Space Telescope / Advan ced Camera for Surveys 
F814W I-band dKoekemoer et alll2007h . 

The observations are being carried out with the 8-meter ESO 
Very Large Telescope using VIMOS, a multi-slit spectrograph 
with a resolution of R ~ 600, resulting in an r.m.s. accuracy on 



the ve locity determinations of ~ 110 km s 1 dLillv et all I2007L 
l2009h 

This paper is based on the first ~ 10,000 spectra, which 
have a rather no n-uniform sampling pattern (see Figure 5 in 
iLillv et al.Ll2009l) over an area of 1 .40 deg 2 and an average sam- 
pling rate of ~ 33%. From this sample, we extracted a "statis- 
tical sample" of 8481 galaxies with high quality redshift deter- 
mination (flag > 1.5) and excluding broad line AGN. Note that 
within our redshift limit of 0.9 adopted in the following analysis, 
the number of broad line AGN is 17, of which 3 are radio loud. 

Absolute magnitud es were computed following the method 
of lllbert et al.l (12005b from the COSMOS photometry (see 
IZucca et all 120091) . which covers a wide range of wavelength 
from the UV to the NIR. This method computes the rest-frame 
absolute magnitudes using the nearest observed band (at the red- 
shift of the galaxy) plus a correction derived from a template fit 
and is the least template dependent possible. 

We used stellar masses (M star )a nd star formation rates 



(SFRs) computed by Bolzonella et aT] d2009l) using Hyperzm ass 
(see details in iBolzonella et all 120091; iPozzetti et all 12007b by 
means of a spectral ener gy distribution fit to the p hotometric 
dat a. The me t hod u ses the Bruzual & Chariot (2003) library and 
the IChabrierl d2003 ) Initial Mass Function. The library provides 
a model for a simple stellar population and its evolution in 220 
age steps. From this database Hyperzmass imposes a set of 10 
exponentially decreasing star formation histories (with e-folding 
times ranging from 0.1 to 30 Gyrs plus a model of continuous 
star formation) in order to fit the photometric data. We asssume 
the values of stellar mass and star formation rate of the best fit 
template. The typical statistical error on the stellar masses is es- 
timated to be ~ 0.20 dex. 

The usually adopted conversions bet ween rad i o lum inosity 
and star formation rate are based on the Salpetei] d 1955b initial 
mass function and for this reason we corrected our logM slar and 
log(SFR) by adding 0.23 and 0.19, respectively, which are the 
median differences of the two quantities computed with the two 
initial mass functions. 

For the radio band we used the VLA-COSMOS survey, a 
large project co nsisting of 23 pointings in an area of 1 .4 x 1.4 
square degrees ( Sc hinnerer et all [2007; Bon di et all [2008). The 
resolution is 1.5 x 1.4 arcsec 2 and the survey, which has a mean 
r.m .s. of ~ 10.5 u J y, dete cted 2501 sources at more than 5cr. 

ICiliegi et al.l (12009b correlated this s ample with the 
COSMOS optical samples (ICapaket alll2007l) . finding 2060 ra- 
dio sources with optical counterparts. We correlate the zCOS- 
MOS "statistical sample" with the VLA-COSMOS sample ob- 
taining a final radio sample of 315 radio galaxies observed spec- 
troscopically out of the 123 1 radio sources with Iab < 22.5. Note 
that our radio data ampling rate is slightly smaller than the av- 
erage zCOSMOS sampling rate because the zCOSMOS survey 
does not cover the entire VLA-COSMOS area. 

Radio luminosities have been computed assuming that the 
radio spectrum is a power law function (defined by S °c v°) with 

a slope a - -0.7. 

We used the density estimtes computed by iKovac et"ail 
(2009) using the ZADE algorithm, which maximizes the statisti- 
cal robustness of the estimators using both the spectroscopic and 
the photometric redshift samples. The algorithm offers the densi- 
ties computed for each galaxy both as a function of the distance 
of the n' h - nearest neighbour and in spheres of fixed comov- 
ing radius. These densities are also computed by weighting the 
galaxies by either their B band luminosities or stellar masses. As 
tracers of the density field, both galaxies of a flux limited and 
volume limited sample have been used. 



Bardelli et al.: Properties and environment of radiogalaxies 



3 



The method has been tested on mock cosmological simula- 
tions, by extracting samples with the same observational limits 
as zCOSMOS. The overdensities were found to be correctly es- 
timated at all redshifts with a slight (~ 20%) underestimate of 
the real value at low densities, corresponding to (1 + 8) < 0.1 
(where 5 — p- < p > / < p >). In the following, we will use the 
densities derived with the nearest neighbor approach. 

I n addition, the zCO SMOS optically-selected group cata- 
log dKnobel et all 120091) of fers a comp l ement ary perspective 
on the role of environment dlovino et all [2009). These groups 
have been extracted using a combination of the friends-of- 
friends method with a Voronoi tesselation density estimator on 
the zCOSMOS spectroscopic catalogue. The method, applied 
to the statistical sample, detected ~ 800 groups, 103 of which 
have more than 5 members. The fraction of galaxies in group 
is ~ 20% with some dependence on redshift. Estimates of the 
group richness have been calculated using the number of galax- 
ies of the volume limited sample, considered as proxies for the 
total richness of groups. 

3. AGN versus Star Forming galaxies in 
environmental studies. 

There are essentially two steps in studying environmental ef- 
fects for radio sources: the division of radio sources according 
to AGN and star forming induced radio emission and the defini- 
tion of coherent control samples. In fact, it is already known that 
radio AGN reside mostly in early type galaxies and therefore 
their overdensity distribution has to be compared only with that 
of radio quiet early type galaxies, otherwise any possible differ- 
ence observed in environme nt could just b e due to the optical 
morphology-density relation dDresslerlll980l) . 

However, a sharp division between AGN and SF galaxies in 
some cases has little physical meaning, because we know that 
both phenomena are commo nplace in a number of galaxies (see 
e.g. ISilverman et all l2009al) and it is therefore difficult to de- 
termine the fraction of the radio flux that originates from AGN 
and/or from the star formation activity. 

For AGN/SF galaxy luminosity fu nction studies (e.g. 
iBardelli et an.l2009tls"molcic et alll2009allbl) . the aim is to count 
the total number of AGN (or star forming galaxies) and for these 
studies it is correct to consider all objects exhibiting AGN (SF) 
activity, regardless of its relative contribution to the radio flux. 
For density studies the situation is slightly different. In this case, 
we ask whether the environment has (or not) any effect in trig- 
gering radio emission from AGN or star formation and it is not 
easy to reach a conclusion in those composite objects, consider- 
ing that the trigger mechanisms are likely very different. 

In other words, if a class of radio sources exhibits an envi- 
ronmental dependence, is this related to its star forming or AGN 
part of the flux? 

In the following, we present a method for selecting sub- 
samples in such a way that the "radio AGN sample" contains 
radio sources for which the AGN emission is at least about one 
order of magnitude higher than that expected from the star for- 
mation. Therefore, our definition of radio AGN is "an object for 
which the radio emission is largely dominated by the emission 
of the central engine". 

With this definition, the optical AGN as defined on the ba- 
sis of emission line ratios, will be included in the star forming 
galaxies sample if they have a low ratio of AGN/SF radio flux. 
This inclusion is suppor ted by a) the inde pendence of emission 
line and radio emission dBestet all 120051) . at least for low lumi- 
nosity AGN as ours, b) the independence of these AGN from the 
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Fig. 1. Upper panel: SFR(radio) versus SFR(SED) for all ra- 
dio galaxies in the initial VLA-COSMOS sample with / < 
22.5 (black circles). In red the 315 objects within the zCOS- 
MOS survey are highlighted. The solid black line is the 
SFR(SED)=SFR(radio) line drawn for reference. Lower Panel: 
Ratio of the SFR as determined in the optical (from SED and 
lines) and the SFR estimated from radio luminosity for a sam- 
ple of "bona fide" star forming galaxies (see text). In red the 
SFRs obtained from the [Oil] line, while in blue from the Ha 
line. Black points refer to SFRs determined from SED. Blue and 
black lines are the least squares fits to line and SED determina- 
tion. 



environment dSilverman et all l2009bh . similar of what has been 
a posteriori found for star forming galaxies and c) the number 
of optical AGN with radio emission, which is negligible with 
respect to that of star forming galaxies. 

3.1. AGN versus star formation radio emission 

There are several recipes for dividing the radio population into 
sources w ith emission f r om an AGN or from star formation (SF) 
(see e.g. ISmolcic et all l2008t iBest et ail I2005I and references 
therein). They are mainy based on determining the optical prop- 
erties of known radio emitting star forming galaxies and defining 
a l imit for th e emis sion of the AGN. In particular, the method 
of IBest et all d2005l) separates the two classes on the basis of 
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SFR [SED] c [M Q yrs" 1 ] 

Fig. 2. SFR(radio) versus the corrected (see text) SFR(optical) 
for our zCOSMOS sample. Red points are star forming galax- 
ies, blue squares are optical AGN (LINERs and Seyfert 2) as 
determined by line diagnostics, while empty squares are radio 
emitting objects in the passive region in the infrared colors- 
specific star formation plane (see Section 3.2). The green line 
is the adopted division between star forming galaxies and AGN. 
Small black points are objects without spectral diagnostic clas- 
sification (due to low quality line measurements). 



the D„(4000)-Li AGHzl^star plane. In practice, this method com- 
pares the ages of stars (estimated with the D„(4000) parame- 
ter) with the current specific star formation rate (estimated with 
L\AGHzlM s tar) and the chosen separation line is 0.225 above the 
D„(4000) value of a 3 Gyr e xponentially decaying star formation 
galaxy track. The method of Smol cic et alj Q2008) is based on the 
fact that the two families define two peaks in a PCA-based rest- 
frame colors combination. Both methods have been calibrated 
with the use of line diagnostics. 

Our method instead has two steps. First, we compare the star 
formation rates estimated from the radio emission and from the 
spectral energy distribution fit dBolzonella et al.L l2009t) and we 
separate AGN from SF radio sources on the basis of this com- 
parison. For the estimate of the star formation rate from the ra- 
dio emission we used the lBelll d2003l) formula. Second, we study 
the optical and infrared properties of the radio galaxies to define 
control samples with similar properties. 

Note that our results do not depend strongly on the specific 
formula adopted in converting L ra </ lo to SFR because our method 
is based on a relative comparison, which considers as AGN all 
objects outside a given radio-optical star formation correlation. 
With the zCOSMOS spectra we could also use the star forma- 
tion rates derived from the lines, but we decided to use the SED 
determinations because: a) at low redshifts (z < 0.5) the emis- 
sion line estimate would use the Ha line (which is affected by 
problems with sky line contamination and strong fringing in the 
spectra in the redshift range [0.3-0.5]), while at higher redshift 
the usable line is [Oil] and b) for weak star formations the line 
measurements are able to provide only upper limits. We note 



that because of point a), the SFR estimate would not be homo- 
geneous being based on two different lines in two different red- 
shift ranges. Moreover, by considering only detected lines, the 
final usable sample would be ~ 60% of that we have with SED 
determinations. For this reason, we will use the spectral features 
only to calibrate the method. 

In Figure Q] (upper panel) we show the SFR(radio)- 
SFR(SED) plane for the entire VLA-COSMOS sample limited 
to / < 22.5) and its zCOSMOS subsample (red points). Only in 
this case, used for illustr atory purpo s e, we used the photomet- 
ric redshifts as derived by Ilb ert et ail (|2009) both for computing 
distances and SEDs. 

It is clear that at high (> 10M G yr" 1 ) values of SFR(SED), 
the two quantities are reasonably correlated, defining a "star 
forming galaxy line", while at low values (SFR(SED)< 2M Q 
yr~') a significant population of radio sources with a radio star 
formation "excess" is present. This excess means that the radio 
emission of this population is not driven by star formation and 
thus is due to AGN. 

However, we need to test the robustness of the SED star 
formation rate determination. By taking advantage of hav- 
ing the zCOSMOS spectra, we can apply standard line diag- 
nostics (Baldwin Phillip s~& Terlevicfl, 1 1 98 lb to classify emis- 
sion line galaxies as star forming, Seyfert 2 or LINERS. This 
procedure wa s applied in an autom ated way to the zCOS- 
MOS sample (j Bongiorn o et all 120091) following the method of 
lLamareille et al] (120091) . At low redshifts (0.15 < z < 0.45) the 
[OIII]/Hy0 versus [NII]/Ho' ratios have been used (red diagnos- 
tics), while at higher redshift (0.5 < z < 0.93) the objects have 
been classified using the line ratios [OIII]/H/3 vs [OII]/HyS (blue 
diagnostics). 

From the 133 radio detected galaxies classified as star form- 
ing, we selected a subsample of 83 "bona fide star forming galax- 
ies", i.e. classified as SF galaxies from diagnostics and with high 
quality line measurements and with emission lines far from sky 
and fringing lines. In FigureQ](lower panel) we plot the ratios of 
the SED (black points) or lines (red and blue points) SFR to the 
radio SFR determination versus the radio SFR. Red points corre- 
spond to SF R obtained from the fOII] line and blue p oints from 
the Ha line (Mousta kas Kennicutt & Tremonti.l l2005). Note that 
SED and emission line determined SFR are consistent with each 
other, but their ratios with the radio SFR shows a significant de- 
crease with increasing SFR(radio). 

In particular, assuming that the values of the radio SFR rep- 
resent the true ones, the SFR derived from SED and emission 
lines overestimate the low and underestimate the high SFR. We 
note that the underestimate of the SFR from the Ha line (red 
p oints) for SFR > 1M yr~' is qualitatively similar to the result 
of ICaputi et al] (120081) . found by comparing this estimation with 
the value derived from the combined UV and IR luminosity (see 
their figure 8). 

We fitted the points in the figure (the "bona fide star 
forming galaxies") with a linear relation finding that < 
/og(SFR[SED]) >= 0.551og(SFR[Radio]) + 0.347, where < 
/og(SFR[SED]) > is the mean value of log(SFR[SED]) at a 
given /og(SFR[Radio]). Then, we computed the "corrected" 
SFR[SED] C by applying the relation 

/og(SFR[SED] c ) = log(SFR[SED])- < log(SFR[SED]) > 

+fc>g(SFR[Radio]) 

. In practice, we force the two estimates of SFR for the 
"bona fide star forming galaxies" to statistically follow the rela- 
tion /cg(SFR[SED] c ) = log(SFR[Radio]) The r.m.s. dispersion 
around the fit is ~ 0.35. 
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Finally, we define the division line between SF galaxies 
and AGN to be the line parallel to the Zog(SFR[SED] c ) = 
log(SFR[Radio]) relation and shifted horizontally by 3 sigma. 
In this way, to be included in our AGN sample, the radio emis- 
sion from the black hole should be about an order of magnitude 
higher than the average emission due to star formation. 

Obviously, this division line between star formation and 
AGN is somehow subjective, because it is likely that there is 
some AGN contribution to the radio flux also in objects below 
the division line. However, this procedure ensures that we define 
a "clean", even if not complete, AGN sample. In particular, at 
higher star formations rates, we cut progressively at higher radio 
luminosities. 

Applying this procedure we define a sample of 97 of AGN 
out of our total sample of 315 radio sources (Figure|2|. 

Among the 133 automatically diagnosed star forming galax- 
ies, only 13 are in the AGN region. After visual checks of the 
spectra, we conclude that all but two could be considered incor- 
rect or at least uncertain classifications because of a low signal- 
to-noise of the lines or contamination by the sky/fringing lines. 
This shows that a large fraction (> 90%) of the objects classi- 
fied as star forming galaxies on the basis of the optical spectral 
diagnostics are indeed recognized as such by our method. 

As expected, the situation for optically classified AGN is 
less clear: among the 14 LINERs and 1 1 Seyfert 2 objects, the 
galaxies showing "star formation excess" are 2 and 3, respec- 
tively. Assuming that for our optical AGN the average contribu- 
tion of the AGN to the l ine luminosity is ~ 50% (see figure 6 of 
ISilverman et all 12009a) and comparing the expected radio star 
formation from the optical-radio star formation correlation with 
the observed radio luminosity, we conclude that, on average, the 
AGN radio luminosity for these objects is a factor ~ 3 - 4 higher 
than that due to star formation. 

Note that although the SFR(optical)-SFR(radio) plot is use- 
ful for distinguishing AGN from star forming galaxies, it can not 
be used to extract a consistent control sample, because the divi- 
sion between galaxy classes is not derived using optical proper- 
ties only. 

3.2. The color-specific star formation plane 

We now discuss the relation betwe en the infrared color 8.0-4.5 
yum obtained by the Spitzer satellite dSanders et al.Ll2007h and the 
optical specific star formation rate (SSFR). It is already known 
that near infrared co lors are a powerfu l method for finding broad 
line AGN (see e.g., lLacv et all 120041) . but here we are more in- 
terested in the behaviour of the "normal" galaxy population. We 
define the infrared color as ratio of the two fluxes S 8.0/S4.5 

In the upper panel of Figure[3]we show the distribution of the 
galaxies of the statistical sample (as defined in Section 2), where 
the star forming galaxies are indicated in red, the Seyfert 2 galax- 
ies in green and the LINER galaxies in blue, all classified with 
spectroscopic diagnostics. Note a number of non-classified ob- 
jects (black points) that reside in the star forming galaxies, where 
the emission lines are expected to be strong. These classification 
failures can be explained by noting that the redshift distribution 
of these objects has two peaks at z ~ 0.5 and z ~ 0.85. In the 
first case, we are at the point in which the Ha line is at the red 
border of the observed spectrum and affected by fringing, and 
the [Oil] line starts to enter from the blue border. At the peak at 
high redshift, the H/3 starts to be affected by fringing. 

The existence of two populations is rather obvious and is 
related to the bimodality found also in the optical by using B-I 
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Fig. 3. Upper panel: infrared color versus specific star forma- 
tion rate for all galaxies in the statistical sample. Red points rep- 
resent the star forming galaxies, blue points the LINERS and 
green points the Seyfert 2 objects classified by line diagnostics. 
Black points represent unclassified galaxies (see text). Lower 
panel: same as the upper panel, showing only the radio galaxies. 
The large circles represent our defined radio AGN. Blue squares 
are the optical AGN (not divided into LINERs and Seyfert 2). 
Vertical and horizontal lines represent the division between the 
three galaxy populations defined in the text. 

or U-B colors (see lCucciati et all l2009) for dividing red, passive 
galaxies from the star forming ones. 

Interestingly, the Seyfert 2 objects occupy the same locus of 
star forming gal axies, while LINERs are at the upper border (see 
also figure 19 of lSmolcic et"aTll2008l) . 
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Passive galaxies are located in a clump in the lower left of the 
plot, while the star forming galaxy region has a vertex nearby the 
passive clump and opens at larger distance. The division value 
has been chosen as the point in SSFR which show a minimum in 
the density distribution of the galaxies. 

In the lower panel of Figure [3] we show the same plot but 
only for radio sources. The large black circles are the radio AGN 
objects defined in the previous subsection. 

This plot shows that the radio AGN sample consists of two 
different populations, one corresponding to the passive galaxies 
(log(S 8.o/5 4 . 5 ) < -0.2 and log(S FR/Mstar) < -11.1) and one 
extending toward the SF galaxies locus. 

Moreover, it also appears as if there is a specific star for- 
mation limit (at Zog(SFR/Mstar) 10.5), such that the radio 

AGN are rare for SSFR higher than this limit. The division at 
log(S FR/Mstar) ~ -10.5 is certainly subjective and may de- 
pends on the characteristics of the survey, but it is a fact that 
the majority of radio AGN are on the left of the line and are not 
spread throughout the plot. This division value has been chosen 
to maximize the number of radio AGN and minimize the con- 
tamination from SF galaxies in the region defined for AGN. 

Of the 75 AGN radio sources with measured infrared color, 
38 are passive galaxies, while the remainig are in the second 
region defined above. We note that this behavio ur is similar to 
that shown in figure 10 of Smol cic~et alj d2008l) using the 8.0- 
4.5 fj.m and 5.8-3.6 /urn color plane. 

In Figure [2] there are a number of radio sources defined 
as passive which are not in the passive AGN sample, because 
they are not radio luminous enough to satisfy our cut. Assuming 
the extreme hypothesis that all these objects are passive AGN, 
this would mean an incompleteness of our sample of ~ 15%. 
Recomputing all quantities with the inclusion of these objects, 
the following results do not change. 

As a general conclusion, we define as radio AGN all objects 
with "radio star formation excess", while "passive AGN" occupy 
the passive (i.e. very low specific star formation) galaxies peak 
of Figure [3] The remaining AGN are defined as "non-passive 
AGN", while all other radio sources are considered to be star 
forming galaxies. 

As outlined above, this choice is justified in the light of the 
analysis of the environmental dependencies. It is worth empha- 
sizing that a more sophisticated method to separate AGN from 
starforming galaxies, when both phenomena are present, is still 
lacking in the literature. Our method leads to a contamination of 
the star forming galaxy sample by AGN of ~ 8% (after having 
applied the limits of Table Q]). We checked that this contamina- 
tion has no effect in the following conclusions. 

Therefore, the control sample will be extracted using the lim- 
its taken from Figure [3] and reported in Table Q] The only excep- 
tion is the additional cut that we applied to the stellar mass in the 
passive AGN sample, which is discussed in Section 4.2. 

3.3. Comparison with other classifications 

To check the consistency between our method for classi- 
fiying AGN and SF galaxies and other methods present in 
the literat ure, we will consider the s pectro scopic line di- 
agnostics dBaldwin Phillips & TerlevichL 1 19811) . the D„(4000)- 
LiAGH z/Mstar plane of iBest et~ai! ( 12005)) and the ISmolcic et al.l 
(2008) method. In Figure |4] we present a plot of the equivalent 
width of the [OIIIJ/H/3 versus [NIIl/Ha (red diagnostics) and the 
[OIII]/Hy6 versus [OII]/H/3 (blue diagnostics) line diagnostics for 
radio SF galaxies as defined by our method. We remind that the 
two diagnostics were applied to two redshift ranges ([0. 15-0.45] 



and [0.5-0.93], respectively). Circles are objects for which all 
emission lines were detected. Here, an upper or lower limit im- 
plies that for each line ratio one line has been detected and for 
the other we have an upper limit. When we have upper limits for 
both lines, no point is shown in the diagram. 

In the upper p anel, the solid line is the theoretical 
iKewlev et alj d2001l) divis ion between SF galaxie s and AGN, 
while the dashed line is me|Ka uffmann et al. (2003) li ne. The re- 
gion b etween the IKewlev et al.l d200ll) and the Kauffman n et al.l 
d2003l) lines defines a region of so-called '"com posite galaxies". 
In the lower panel, the division proposed by Lamareille et al.l 
(2004) is plotted as solid line and dashed lines represent the 
+/ - 0.15 dex uncertaintes in the classification. 

Of the 100 SF galaxies within the redshift range useful for 
the red diagnostics, 8 1 have spectra with all four lines detected, 
while for the blue diagnostics the ratio is 63/91. As expected 
from Figure |2l a not negligible number (14) of the optical AGN 
(i.e. those defined by line ratios) are in the SF galaxy sample. 

For the red diagnostics, 18 objects fall in the "composite 
galaxy region": by applying the same procedure defined in Sect. 
3.1 to estimate the relative contributions of AGN and star for- 
mation rate to the radio flux, for these objects we obtained on 
average ~ 2.5. For the non-passive and passive AGN samples, 
the diagnostics are inconclusive, because of the lack of detected 
lines for estimating the ratios. The blue and red diagnostics give 
a different ratio between the number of detected AGN and star 
forming galaxies. This is related to the different efficiency in de- 
tectin g AGN of the two diagnostics, as noted in Bong iorno et al.l 
(2009) (see their Figure 3 and 4). In particular, the blue diag- 
nostics classify a large number of X-ray emitting AGN as star 
formimg galaxies, reinforcing therefore the need of new AGN- 
star forming division methods. 

Another method for differentiating between AGN and SF 
galaxies through the use of the D„(4000)-Li ACHz/M slar plane 
dBest et aU 120051) . where the age of the galaxy stellar popula- 
tion (estimated by D„(4000)) is compared with the specific star 
formation rate derived from the r adio luminosity. The division 
line adopted bv lBest et all d2005l) is 0.225 above the D„(4000) 
value of the track defined by a 3 Gyr exponentially decaying star 
formation. In Figure [5] we show our samples in this plane. The 
number of AGN below the division lin e is 5 (out of 104) , while 
~ 50% of our SF galaxies are above the lBest etal] (f2005) divid- 
ing line and they would be classified as AGN on the basis of this 
criterion. These galaxies are likely in a transition region, where 
our AGN induced emission objects are mixed with the SF ones. 

The ISmolcic et al.l d2008l) PCA based classification is an- 
other, complementary way for classifyi ng AGN and star for m- 
ing objects. Of the 129 AGN defined in ISmolcic etalJ d2008l) in 
common with the zCOSMOS sample, 66 are star forming ex- 
cess galaxies in our class ification scheme (and therefore AGN), 
while 9 out of 97 of the Smol cic et alJ (2008) star forming ra- 
dio galaxies are in our AGN locus. The fraction of AGN that 
have log(SSFR) < -10.5 is ~ 60%. A similar percentage of 
st ar forming excess wa s found in the AGN sample as defined 
in iBardelli et al.l (2009), where a spectrophotometric division is 
adopted. The appar ently large difference between our and the 
Smolcic et al. (2008) method is largely due to the fact that opti- 
cal A GN and composite galaxies are included in the AGN sam- 
ple in Smolcic et al. (2008) and not in our scheme. 

In fact, by applying our diagnostics to the ISmolcic et al.l 

(2008) sample, we confirm that the contamination of their AGN 
and star forming sample is ~ 17% and ~ 10%, respectively con- 
sistent with what is estimated in their paper. 
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Table 1. Sample definition: numbers in the third and fourth columns are with and without the "evolving" magnitude limit 



sample 


limits 


Number in radio sample 


Number in Control Sample 


"passive AGN" 


log(SFR/Mstar)[SED]<-ll.l & log(5 8 .o/5 4 .5) <-0.2 


38/38 


524/530 


"non-passive AGN" 


-11.1 <log(SFR/Mstar)[SED] <- 10.5 


31/37 


496/581 


star forming galaxies 


log(SFR/Mstar) [SED] >- 1 0.5 


154/186 


1672/3154 



A general conclusion of this comparison is that our approach, 
which entirely relies on the comparison of observed radio power 
with the expected power due to the star formation, is more con- 
servative (and therefore not complete) in selecting radio AGN 
than other methods. We remind that our method was constructed 
to classify as AGN those objects for which the radio flux is dom- 
inated by the central engine. 

4. Differences between radio loud and control 
samples 

4. 1 . Spectral differences between radio loud and control 
samples 

The most basic approach is to construct composite spectra of the 
samples. 

In Figures [6] and [7] we show the composite spectra for the 
passive, non-passive radio AGN and star forming galaxies and 
their control samples. Moreover, we plot the difference between 
the radio and control sample composite spectra. 

For passive AGN, there is no significant difference between 
the composite spectra of the radio loud and control sample, im- 
plying that there is no obvious AGN signature in their optical 
spectra. There is only an indication for radio spectra to be redder 
than the control sample. 

The spectra of both the non-passive AGN and control sam- 
ple objects exhibit significant [Oil] emission lines, as well as 
absorption lines (like H+K Ca and Mg 5175 A) and relatively 
weak [OIII] and Ha lines. 

In both spectra Hfi is in absorption only. All these spectral 
features imply that these galaxies are of relative early type. All 
galaxies from the non-passive AG N and control sample are in- 
deed classified bv lZucca et al.l d2009f) as spectrophotometric type 
1 and 2 (roughly corresponding to Sa-Sb). However, non-passive 
AGN tend to be of even earlier spectrophotometric types than the 
control sample. The radio sample has stronger [Oil] equivalent 
widths (15 A versus 7 A) with respect to the control sample. 

Considering the star forming galaxies, the difference be- 
tween the two composite spectra is negative in [Oil] and posi- 
tive in Ha, implying that in radio emitting star forming galaxies 
a larger amount of line extinction (from dust) is present. This 
is confirmed by the fact that the continuum of radio emitting 
galaxies is redder than the control sam ple. This is related to the 
differences in colors already found by Bar delli et al.l d2009f) for 
VVDS radio galaxies and present also in our VLA-zCOSMOS 
sample. 

Additional information about the spectral features can be 
derived from the equivalent width distribution of the lines. We 
checked this point by examining the distribution of the equiv- 
alent width of the [Oil] 3727 A line and the 4000A break (as 
measured in Mign oli et all |2009). There are no significant dif- 
ferences between the distributions of these quantities for passive 
AGN and their control sample. 

Viceversa, the radio detected star forming objects have a 
4000 A break slightly shifted toward higher values than the cor- 



responding control sample. Given that the stellar mass distribu- 
tion of the radio detected star forming galaxies and control sam- 
ple is different and that our data show a stellar mass- 4000 A 
break relation, this difference is likely due to the different stellar 
masses sampled (see Section 4.2). In fact, the differences disap- 
pear considering objects with stellar mass > 2 x 10 10 M , for 
which the mass distributions of the radio star forming galaxies 
and their control sample are identical (see Section 4.2). 

More subtle is the difference for the non-passive AGN sam- 
ple (Figure[8]l. In this case the radio detected objects have higher 
values for the 4000A break with respect to the control sample 
and lower values for the equivalent width of the [Oil] 3727 A 
line (the probability for the KS test that the two distributions 
come from the same distribution is ~ 0.005). 

The last question is whether the control sample of the non- 
passive AGN shows some sign of AGN activity in their emission 
lines. In Figure|9j we show the diagnostic diagrams for this sam- 
ple. The blue and magenta circles correspond to the values ob- 
tained from the composite spectra of non-passive AGN and their 
control sample, respectively. Considering the galaxies with both 
all lines detected or only one line detected in the considered ra- 
tio, ~ 50% of this sample consists of optical AGN or upper limits 
consistent with that of AGN. Note that this fraction is higher in 
the red, low redshift diagram. These fractions are broadly con- 
sistent with those of the radio sample, although the number of 
objects with all detected lines or meaningful upper limits is only 
20 for the radio detected non-passive AGN sample. 

We note that the points corresponding to the composite spec- 
tra are not in the middle of the distribution of the ratios of single 
galaxies. This is because the composite spectra are obtained also 
considering objects not plotted in Figure|9]( which are ~ 48% and 
27% of the samples in the blue and red diagram, respectively) be- 
cause of the lack of detected lines. Therefore, composite spectra 
are not expected to have the average ratios of the plotted points. 

4.2. Stellar Mass distributions 

Before studying the effect of the different environments on the 
radio galaxies, it is necessary to verify the consistency between 
the stellar masses of the various samples of radio galaxies de- 
fined in the previous sections with those of the ir respective con- 
trol s a mples. Given the mass- density relation (Scodeggio et all 
2009; iBolzonella et all 120091) . a difference in the stellar mass 
distribution may induce a difference in the observed overdensity 
distributions, which would not be related to AGN or star forma- 
tion mechanism. 

To consider a homogeneous population at all redshifts, we 
selected galaxies and radio sources within an optical luminosity 
complete, volume limited sample at z < 0.9. To take into accoun t 
the observed optical luminosity evolution ( IZucca et al.L [2009), 
the absolute magnitude limit was defined to be Mb = -19.8 - z. 
This limit in absolute magnitude corresponds to the zCOSMOS 
limit in apparent magnitude (/ < 22.5) at z — 0.9. In Figure [TUl 
we show the observed mass distibution of the passive AGN and 
their control sample. 
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Fig. 4. Red and blue diagnostics diagram for our star forming 
galaxies with radio detection. Circles are objects with all lines 
detected. For the upper or lower limits, we plot only objects with 
only one upper limit determination for each line ratio. 



For the control sample, since we have a volume limited sam- 
ple, we plot the stellar mass histogram. Since the mass complete- 
ness limit is changing with redshi ft (from ~ 10 9 M G at z ~ 0.3 to 
~ 5 1O 1O M at z ~ 0.9)(see IPozzetti et all l2009b . our observed 
stellar mass distribution shows an increasing incompleteness at 
decreasing stellar mass. Assuming that the incompleteness is the 
same for the control sample and the radio sample, only the rela- 
tive comparison between the control and radio sample is relevant 
here. The radio sample, which also has the limit in radio flux in 
addition to the limit in apparent magnitude, needs to be corrected 
for V max . 




to 12 IQW 
l0 g( L 1.4GHz/ Mstar ) 

Fig. 5. Dn(4000) versus Li 4G// Z /Mstar plane: Triangles are our 
star forming galaxies and circles our AGN sample. Open squares 
indicate the passive AGN. 



The first result is that all radio emitting objects have 
log(Mstar) > 10.7, while the control sample has a broader mass 
distribution. For this reason, in the following analysis we limit 
the control sample to log(Mstar) > 10.7 . 

We note that the ratio of the stellar mass distribution of ra- 
dio emitting passive AGN to its control sample is rather con- 
st ant at ~ 20%. This is a pparently inconsistent with the results 
of ISmolcic et a l. (2009a), for which the ratio of radio AGN to 
the control sample is rapidly increasing with stellar mass. The 
difference could be understood considering that we study only a 
specific class of AGN and apply a significantly high cut in ab- 
solute magnitude to create a volum e limited sample at the depth 
of zCOSMOS. On the other hand. ISmolcic et alJ d2009a) con- 
sidered volume limited samples in terms of the radio limit (i.e. 
fixing a limit in the radio luminosity) and corrected for V max the 
optical incompleteness. Their apparent magnitude (and conse- 
quentely absolute magnitude range) limit is deeper than that of 
zCOSMOS survey. In fact, considering all AGN (passive+blue) 
and av oiding the absolute ma gnitude limit, we derive similar ra- 
tios as ISmolcic et al.l (l2009al) . 

For non-passive AGN the ratio of the two stellar mass distri- 
butions (radio detected and control sample) is ~ 5% for masses 
< 10 11 M (see upper panel of Figure[TT), while at higher masses 
the ratio is > 25% We note that in this case we recovered the 
mass dependence of the AGN activity. 

For star forming galaxies, the lower panel of Figure [TT] in- 
dicates that the the two stellar masses distributions are the same 
for masses higher than 2 x 10 10 M B and differ for low er masses, 
consistent with a similar plot in Bardell fet al.1 (120091) . 

We interpret this trend as implying that the stellar mass dis- 
tribution of radio detected star forming galaxies is the convolu- 
tion of the stellar mass function of these objects with the star 
formation r ate distribution at fixe d stellar mass or absolute mag- 
nitude (see iBardelli et all 120091) . The mean star formation rate 
increases with mass, but the distribution has significant tails. At 
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Fig. 6. Composite spectra for the radio emitting AGN (in black, 
upper spectrum) and the control sample (in red, middle spec- 
trum). The upper panel refers to passive AGN and the lower 
panel to non-passive AGN. The spectrum of radio sources is 
shifted for clarity by 0.9 in the y-axis. At the bottom of the fig- 
ures, we plot the difference between radio emitting objects and 
control sample (lower spectrum). 



lower stellar massses an increasing number of objects with star 
formation rate in the low tail is below the flux limit of the radio 
survey. More specifically, the lower limit for the star formation 
rate is set to fc>g(SFR/Mstar) ~ -10.5 from Figure[3] which cor- 
responds for log(Mstar) = 10.3 (the "completeness limit") to 
~ 0.6M Q yr~\ similar to the minimum observed star formation 
rate. For lower masses the limit is lower than the minimum ob- 
served star formation rate. 
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Fig. 7. Same as Figure [6] but for the radio emitting star forming 
sample and its control sample. 
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Fig. 8. Equivalent width and Dn(4000) distributions for the radio 
and control sample for the non-passive AGN. The dashed red 
line is for the control sample and the black solid line is for the 
radio sample. 



For the density analysis we limit the control sample to 
log(Mstar) = 9.7 for the star forming galaxies and to 
log(Mstar) - 10.4 for non-passive AGN. These limits corre- 
spond to the minimum mass detected in the radio samples. 
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Fig. 9. Line diagnostics for the control sample of the non-passive 
AGN sample. The blue and magenta points correspond to the 
values obtained from the composite spectra of non-passive AGN 
and control sample respectively. 



5. Environmental effects 

5.1. Method 

To consistently define our radio and control sample we adopted 
the limits reported in Table Q] We computed the density distri- 
butions by considering all objects with z < 0.9 and brighter than 
an "evolving" absolute magnitde limit Mb < -19.8 - z, where z 
is the galaxy redshift. This evolving limit was introduced in or- 
der to take roughly into account the luminosity evolution of the 
galaxies as determined fr om the luminosity function estimate by 
IZucca etail J2009, 2006) and therefore to approximately con- 




Stellar Mass [M ] 

Fig. 10. Stellar mass distribution of the passive AGN (black con- 
tinous line) compared with its control sample (red dashed line). 
The vertical dashed line indicates the adopted stellar mass cut. 
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Fig. 11. Upper panel: Stellar mass distribution of the non-passive 
AGN (black continous line) compared with its control sample 
(red dashed line). Lower Panel: same but for star forming galax- 
ies. 



sider the same population at all redshifts. Moreover, in the case 
of the passive AGN sample, we also applied a cut in stellar mass 
at log(Mstar) > 10.7, for the non-passive at log(Mstar) > 10.4 
and for the star forming galaxies at log(Mstar) > 9.7 (see 
Section 4.2). 
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5.2. Redshift distribution 

The most basic approach to assess a possible dependence of ra- 
dio galaxies on the environment is to compare their redshift his- 
togram with that of the control sample. In Figure [12] (upper and 
lower panels) the normalized histograms of the radio sources 
(solid black line) and of the control sample (dashed red line) 
are shown. 

In this figure, we cut the radio sample at logL.R fl </, (W Hz _1 )> 
23 to obtain a volume limited sample also in terms of radio lumi- 
nosity. Variations in the ratio of the number of radio galaxies to 
those in the control sample found in correspondence to peaks or 
valleys of the galaxy distribution would indicate a different en- 
vironment distribution between the two populations. Smoother 
variations as a function of redshift would imply an evolution of 
sources. 

In the upper panel we consider all AGN (non-passive + pas- 
sive), i.e. those radio sources with star formation rate excess and, 
for the control sample, all galaxies with log{S S FR)(optical) < 
-10.5. In the middle panel, we show the ratio of the two his- 
tograms in order to visualize better any redshift dependence. 
Errors are presented at the 1 <x confidence level. Black points 
correspond to the whole AGN sample, while red points refer to 
passive AGN. For clarity the red points have been shifted by 
+0.03 along the x axis. 

In the lower panel we show the redshift distribution of star 
forming galaxies. 

The number ratios of both AGN and SF radio samples to 
the respective control samples are rather constant. This implies 
that the radio samples follow the same large scale distribution 
as the control samples and therefore no large segregation with 
the density is present at the largest scales. For the two classes of 
AGN taken separately, the ratio values are around 0.07 for the 
non-passive AGN sample and 0.06 for the passive AGN sample. 

The redshift histograms explore densities which are on the 
largest scale possible for our redshift survey, ranging from ~ 12 
Mpc at z ~ 0.3 to ~ 21 Mpc at ~ 0.5. In the following, we 
explore smaller scales using the 5' /! nearest neighbor. 



5.3. The nearest neighbour determined overdensity 
distribution 

We used the de nsities estimated fro m the distance of the 5'' 1 near- 
est neighbour dKovac et al.L 12009) and the volume limited sam- 
ple as tracer of the density field. In this case, the tracers are the 
galaxies with absolute B band magnitude brighter than — 19.5 — z 
up to redshift 0.7, while for redshifts between 0.7 and 1, the 
tracers are galaxies with Mb < -20.5 - z. This corresponds to 
overdensities computed on average scales of ~ 2.4 Mpc, thus 
not much larger than the typical extent of a clus ter. We also 
checke d the results for the other tracers described in Kova c et ail 
(2009), but the results did not change. Among the various avail- 
able weighting schemes, we used the densities weighted by the 
stellar masses of the galaxies, but the results did not change ei- 
ther when using the no-weight (i.e. number weighted) estimator. 

In Figure Q~3] we compare the overdensity distribution of the 
three radio samples (black solid histograms) with those of their 
respective control samples. In the upper left panel, we present the 
total AGN sample. We recall that in this case the radio sample 
comprises all radio sources with star formation excess and the 
control sample comprises all objects with log(S S FR)(optical) < 
-10.5. The two distributions are different at about 2cr level on 
the basis of the K-S test. 
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Fig. 12. Upper and lower panels: fractional redshift distribu- 
tions of our samples. Red dashed histogram: control sample. 
Black solid histogram: radio sample. Histograms have been nor- 
malized to the number of objects. Upper Panel: AGN sample 
(passive+non-passive). In the middle panel, we show the ratio 
between the number of radio AGN and the control sample. Red 
points represent the fraction of radio loud passive AGN (shifted 
for clarity by 0.03 in the x axis). Lower panel: Star forming 
galaxies redshift histograms and relative ratios (black points). 



A much stronger statistical signal of difference (as evaluated 
by the -KS test probability) is obtained by considering only the 
passive AGN (upper right panel). Here, we note a lack of radio 
sources at low overdensities (< 3) and an excess at high overden- 
sities (> 10). In fact, the signal present in the whole AGN sample 
is coming entirely from passive AGN. If we select only the non- 
passive AGN sample, the KS test does not find any significant 
difference (lower left panel). In Figure[l4]we plot the same data 
but now as ratio of the number of radio emitting passive AGN to 
the control sample. As can be seen, the ratio steadily increases 
with the density from 0.02 to 0.2. 

There is the possibility that this behaviour is caused by the 
fact that the observed ratio of the stellar mass functions of ra- 
dio to the control sample shows a strong variation between 0.02 
(in the bin at lower mass in Figure ITOb to ~ 0.2 (in the high- 
est bin). Note that in Figure [10] we showed the Vmax corrected 
mass distributions. In other words, spurious density differences 
caused by the stellar mass-density relation could be induced both 
by real differences in stellar mass fucntions between radio and 
control samples, but also by the fact that the observed mass dis- 
tributions are different (although the mass functions are similar) 
merely because of the flux limit. 

To check this fact, we recomputed the control sample over- 
density histogram of the upper right panel of Figure [13] by 
weighting all masses by the ratio of the observed mass distri- 
bution of the radio to that of the control sample. The results 
changed by a negligible amount. We checked a possible depen- 
dence with the redshift of the difference between the environ- 
ments of radio and control sample. The differences between ra- 
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dio and control samples remain significant although the density 
distributions of both samples have higher density tails at lower 
redshift, consistent with the evolution of galaxies and of the 
overdensities. We conclude that no redshift evolution is present 
in the density-radio emission relation. 

Another apparently significant difference in environment 
distribution between radio and control sample is coming from 
the star forming galaxies. The radio emitting star forming galax- 
ies reside in systematically higher densities than the control sam- 
ple. However, we note that also for star forming galaxies the stel- 
lar mass functions for the radio and control sample differ (lower 
panel of Figure fTTb. 

We repeated the previous procedure to correct the control 
sample for the stellar mass distributions ratios and we find that 
the significance of the difference disappears. To do this we com- 
puted the x 1 probability of the two distribuions: we can not use 
the KS test because of the binned mass ratios assumed for the 
correction: the probability that the two distributions are the same 
increases from ~ 0.05 to ~ 0.90. 

Therefore, we conclude that most of the difference between 
the overdensity distributions of radio and control star forming 
galaxy samples is caused by the different stellar masses distribu- 
tion sampled by radio and control samples (higher masses reside 
in higher densities) and not by the physical role of the environ- 
ment on the radio emission. 

All results showed in this paragraph are stable also using the 
10" 1 and 20"' nearest neighbor, which correspond to scales of 3.5 
and 5 Mpc, respectively. 

5.4. Groups of galaxies 

In Figure [15] we show the ratio of the number of radio pas- 
sive AGN to the control sample as a function of group richness. 
The fraction of passive galaxies which are AGN in the isolated 
galaxies sample is 0.052 (16/307), while the fraction in groups is 
0.104 (14/135), i.e. the difference is not statistically significant. 
It is however worth noting that 2/3 of the control sample galaxies 
are isolated and 1/3 are in group, while the radio passive AGN 
are nearly equally distributed in the two classes. 

However, if we limit the group sample at richness> 4 (for 
which the contamination by false groups is small) the fraction 
of radio passive AGN becomes 0.138 and the significance of 
the difference between the fraction of passive AGN among iso- 
lated galaxies and in groups to 2.15cr. In fact, the overall distri- 
butions of radio and control samples has a significance derived 
from the K-S test of 8% to be drawn from the same distribution. 
Therefore, there is a small signal coming from an increase of the 
fraction of radio objects with richness. 



6. Properties of the passive AGN 

In the previous Section we showed evidence that passive AGN 
have an increasing probability to be a radio source in denser en- 
vironments. Here, we explore if other properties of single pas- 
sive AGN also change as a function of local density. In practice, 
we divided the sample of passive AGN in two bins at a value of 
1+5=10, and compared a number of quantities for the two sub 
samples. First of all we computed the two luminosity functions, 
We note that in our case it is only possible to consider differ- 
ences in shapes, as the relative normalizations are influenced by 
the volume occupied by the various overdensities within the sur- 
vey volume, and this is difficult to control. For this reason, we 
compared only the shapes which resulted to be indistinguishable. 




Fig. 14. Ratio between the radio and control sample of passive 
AGN as a function of the density. 




1 10 

log(Richness) 



100 



Fig. 15. Fractional distribution of richness of the groups where 
radio passive AGN (solid line) and galaxies of the control sample 
(red dashed histogram) reside. Circles with error bars represent 
the ratio between the numbers of objects of these two samples. 



This is not unexpected because already ILedlow & Owed fl996) 
found that the radio and field luminosity functions have similar 
shape. 

In Figure [16] (upper panel) the radio power-stellar mass ra- 
tio is plotted. Dashed red lines correspond to the high density 
sample and the solid blue line is the lo w density one. Each ra tio 
has been weigthedby its V max following Bardelli et al. (2009). It 
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Fig. 13. Black histograms: overdensity distribution of the radio sources, Red dashed histograms: control samples. Upper left panel: 
total AGN sample. Upper right panel: Passive AGN sample. Lower left panel: non-passive AGN. Lower right panel: Star forming 
galaxies. The blue histogram corresponds to the corrected control sample. 



appears that, at a given stellar mass, the radio luminosity of the 
AGN is significantly higher in dense environments. 

In the lower panel of Figure[l6] we plot the radio luminosity 
versus the stellar mass. We note that the difference between the 
two panels is that in the upper one the radio luminosity-stellar 
mass ratio is corrected by V max , while in the lower panel we plot 
the observed points. It seems that there is a transition luminos- 
ity ( at logL ra£ / lo (W Hz ')= 23.5 ), below which the objects are 
about equally distributed in low and high density regions, while 
above this luminosity ~ 80% of radio sources reside in high den- 
sity regions. There is also an indication that below the luminosity 



threshold, a correlation between mass and luminosity is present 
(confirmed by a Spearman rank correlation test), which disap- 
pears at higher luminosities. 

In Figure [17] we show the black hole mass distribution 
for the passive AGN (black) and the control sample (red). 
We estimated the blac k hole mass using the prescription of 
iMarconi & Huntl d2003l) 
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Fig. 16. Upper Panel: Distribution of Ls a j, /Mstar for passive ra- 
dio AGN in two bins of overdensity (red dashed line for high and 
blue continous line for low overdensity). Lower Panel: Radio lu- 
minosity versus Stellar Mass for the same samples. Red triangles 
represent high densities and blue points the low densities. 



where Lk is the K-band luminosity of the galaxy. We show 
the re sults for the local Mbh~ L% relation, because Peng et ail 
(2006) claim that th is relation does not cha nge with redshift. We 
also corrected the Marconi & Huntl (120031) formula as done by 
Smo lcic et al.l (l2009a[) for the evolution i n luminosity of a pas- 
sive population (see lHopkins et afl,l2006l) but the results did not 
change. 

There is no significant difference in the mass distribution be- 
tween high and low density AGN. The only difference is between 
the high and low density control sample (with a KS probability 




Black Hole Mass ( M ) 

Fig. 17. Black solid histogram: high density (1 +6 > 10) black 
hole masses distribution for passive AGN. Black dashed his- 
togram: same distribution but for objects with (1 + S < 10). In 
red the distribution of masses for the control sample, shifted for 
clarity along the vertical axis by 0.3. 



for the two distribution to be equal of ~ 0.005). If we plot the 
radio luminosity versus the black hole mass (not presented here) 
we obtain a plot similar to the lower panel of Figure [16] with the 
stellar mass substituted by the black hole mass. This is reason- 
able considering the relation between stellar mass and K-band 
luminosity. 

Our general conclusion is that the main difference between 
passive AGN in low and high environments is a difference in 
radio power and that this is due, given the similarity in the black 
hole masses, to the host galaxies or environment structure. 

7. Discussion and conclusion 

Studying the environment of radio galaxies requires a division 
of the radio sources on the basis of their emission mechanism, 
because AGN and star forming galaxies may be influenced dif- 
ferently by the local density. To do this, we compared the star 
formation rate estimated from the spectral energy distribution 
with the one derived from the radio luminosity. We found a clear 
correlation between these two quantities for the star forming 
galaxies: we defined as AGN those objects for which the radio 
star formation rate is significantly higher than the optical one. In 
other words, these objects have a radio luminosity significantly 
higher than that predicted by their star formation. We assume 
a stringent cut at which the AGN radio emission is more than 
an order of magnitude higher than the star formation one in the 
same galaxy. In this way we excluded from our radio-based AGN 
sample those AGN defined on the basis of their optical emission 
lines, such as Seyfert 2 and LINERs, for which the ratio of AGN 
to star forming radio emission is lower than the adopted thresh- 
old. 

Since AGN are hosted primarly by early type optical galax- 
ies and star forming objects are basically spiral galaxies, we ex- 
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pect to recover at least the optical morphology-density relation 
dDressler[ ll980), not related to the radio emission. Therefore, to 
investigate whether there is an environmental effect on the ra- 
dio emission, we defined control samples to be compared with 
the AGN and star forming radio samples. To achieve this, we 
studied the positions of radio sources and zCOSMOS galax- 
ies in an infrared color-specific star formation plane. We de- 
fined three populations: passive AGN, non-passive AGN and star 
forming galaxies. The passive and non-passive division is rem- 
in escent of the lo w and high-excitacion dichotomy described 
in Smolcic (2009) and probably corresponds to two different 
processes of accretion (hot and cold, respectively). In fact, the 
two proce sses are expected to depend differ ently on the envi- 
ronment (Har dcastle Evans & Crostonl 120071) . The radio lumi- 
nosity functions of the passive and non-passive AGN are sim- 
ilar in shape but not in luminosity range. Non-passive AGN 
stop at logL ra </ IO (W Hz *)= 23.4, while the passive AGN reach 
logL radio (W Hz" 1 )= 25.2. 

Finally, o ur radio objects belongs almost all to the FR-I class 
(as defined bv lLedlow & Owerll 19961) and the sample lack of the 
broad line AGN class. Therefore we focussed mainly to objects 
triggered by gas accretion and did not considered other mecha- 
nism as e.g. mergers. 

Passive AGN: This population of AGN is hosted by red 
galaxies. Using the morphological classificatio n derived from 
HST/ACS images throug h the ZEST algorithm (IScarlata et al.L 
2007t iTasca et all [2009), we find that ~ 85% of both passive 
AGN and relative control sample are classified as early type (i.e. 
bulge dominated). The spectra of these galaxies with and without 
radio emission are indistinguishable and the ratio of their stellar 
mass functions is rather constant at ~ 0.2 for masses higher than 
5 x 10 10 M . This value can be interpreted as corresponding to 
the average fraction of time during which the central engine is 
switched on. Note that our ratio is derived comparing the mass 
functions: if we had used the observed number of galaxies, we 
woul d have found values more similar to those reported in Table 
4 of lShabala et all (|2008) ( ~ 0.5 at high masses, ~ 0.15 at inter- 
mediate masses and ~ 0.01 at low masses). 

This population exhibits a strong dependence on the envi- 
ronment: the fraction of objects with radio emission increases 
from ~ 0.02 for underdense regions (1 + 6 < 3) to ~ 0.15 for 
overdense regions (1 + 6 > 10). Similarly, such trend is also 
found (but at a low significance level) as a function of richness 
of galaxy groups. 

Studying the radio luminosity-stellar mass plot, we found 
that for luminosities lower than log P(W Hz~')=23.5 there is a 
correlation between the two quantities, that disappears at higher 
powers. The correlation between stellar mass and radio lumi- 
nosity is not unexpected considering the relation between the 
accretion on to a black hole in the "radio mode " to the black 
hole m ass and hot g as fraction (see equa tion 10 of lCroton et all 
2006). As stated by ICroton et all (120061) . the hot gas fraction is 
approximately constant for galaxies with v,,,> > 150 km s _1 and 
therefore the black hole luminosity scales with the black hole 
mass, which correlates with the stellar mass. Interestingly, half 
of the objects in the low radio luminosity regime reside in low 
densities, while the others are in overdense regions. 

The correlation disappears for radio luminosities higher than 
log Lr ,/, (W Hz~')=23.5. Almost all of these radio luminous ob- 
jects reside in overdense region. On the other hand, no difference 
is present in the black hole mass between high and low densities 
and therefore it seems that in richer environments the emission 
mechanism is more efficient than in low density regions. 



Therefore, we can tentatively say that at low density the 
emission is determined by the host galaxy stellar mass and at 
high density by the host structure (group/cluster) in which the 
galaxy resides. The higher luminosities in denser environments 
could be caused by a greater fuel supply (due to the central cool- 
ing of the gas) and/or to the fact that also outside the galaxy there 
is a significant gas density to confine the radio jet. On the other 
hand, we already know that ellipticals at the center of galaxy 
groups are richer in h ot gas than isolated and non-dominant ones 
dHelsdon et alll200"ll) and this implies that a mechanism driving 
gas within dominant galaxies is or has been active. 

This can also justify a longer time spent by the central en- 
gine in the "switched on" status, as suggested by the increasing 
ratio of passive AGN to control sample as a function of density. 
It would be interesting to investigate whether the high density 
AGN have luminosities that correlate with the central cooling 
flow, following the r elation betw e en rad io luminosity and accre- 
tion mass found by Mit tal et al. (120081). Note that our findings 
are co nsistent with those of both Best! ( 2004b . iMandelbaum et alj 
(120091) . 

Finally, the difference in overdensity between radio-loud and 
control samples remains constant across the studied redshift 
range, implying that no significant evolution in this phenomenon 
is present. 

Non-passive AGN: This is a less well-defined family. In prin- 
ciple, in this class one can find not only AGN but also objects 
such as post-starburst galaxies. In this case, radio emission dis- 
appears on a time scale of 10 8 years, which could be longer than 
the disappearing time of the star formation signatures in optical. 
However, we expect that these objects are rare because of the 
implied short time sca les. In fact, on l y two of the radio detected 
post-starbursts of the Verg ani et al.l (|2009) sample and one of 
the upper limits lie within our AGN region. This means that 
post-starburst galaxies contaminate our sample by 4%. The non- 
passive AGN population occupies a small region of the infrared 
color-specific star formation rate plot, which is approximately 
equivalent to the green valley defined in the optical bands. These 
galaxies are objects with spectrophotometric type earlier than 
Sa-Sb and applying the standard diagnostics based on the emis- 
sion line ratios we found that 50% are classified as optical AGN 
and the remaining are star forming galaxies. The ZEST mor- 
phologies of non-passive AGN are equally divided between early 
types and spirals (with a ~ 8% of irregulars), while the control 
sample has percentages of ~ 60% and 40%, respectively. The 
radio emitting objects of this class follow the same environment 
distribution of the corresponding control sample galaxies and do 
not show the luminosity-density behaviour. 

Star Forming galaxies The composite spectra of radio emit- 
ting star forming galaxies is redder than the control sample with 
a smaller [Oil] and higher Ha equivalent width. 

Considering that the radio detected star forming galaxies are 
objects where the star formation rate is on average higher than 
the undetected ones, this implies that the dust content is higher. 
This increase of dust ext inction with increasin g star formation is 
already known (see e.g. iPannella et all 120091) and could justify 
the relation between the emission lines and the radio SFR esti- 
mators shown in the lower panel of FigureQ] The ZEST morpho- 
logical classification shows that the radio detected objects tend 
to be of earlier type than the control sample, as found for the 
spectrophotometric types. 

The radio emitting star forming galaxies show an appar- 
ent difference in density distribution with respect to the con- 
trol sample, but we have shown that this is due to the differ- 
ent stellar mass distributions of the two samples. Correcting the 
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control sample with the fraction of radio emitting objects, no 
differences in environment are present between radio and con- 
trol sample galaxies. T his is consistent with the conclusion of 
iKauffmann et al.l ((2004) that the star formation does not depend 
on the density for scales larger that the Megaparsec. 
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